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Abstract 

We analyse the predictions for the inclusive branching ratio for B — > X s ~j in a 
class of string-inspired SUSY models with non-universal soft-breaking A-terms. These 
models are particularly interesting since the non-universality of the A-terms plays an 
important role in providing new significant contributions to CP violation effects in kaon 
physics while respecting the severe bounds on electric dipole moments. We show that 
b — > S7 do not severely constrain the non-universality of these models. In particular, 
at low tan/3 (tan/3 ~ 2), our predictions are close to the universal case. For large 
tan/3 (tan/3 > 15), the effect of non-universality is enhanced and stronger constraints 
hold. We find that the parameter regions which are important for generating sizeable 
contribution to e'/e , of order 2 x 10~ 3 , are not excluded. 
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1 Introduction 



Recently there has been a growing interest concerning supersymmetric models with 
non-universal soft-breaking terms The main theoretical reason is that some 

string-inspired models naturally favour SUSY models with non-universality in the 
soft-breaking sector 10-0. Within this class of string-inspired models particularly in- 
teresting are that ones with non-universal A-terms HQ. This is mainly due to the 
relevant role they can play in solving the SUSY CP problem, especially in the light of 
the recent experimental results on the direct CP violation parameter e'/e . 

The new measurements of e'/ e at KTeV 0] and NA48 [llj lead to a world average 
of Re e'/e =(21.4±4.0) x 1(T 4 p[. This result is higher than the Standard Model (SM) 
predictions |13|, opening the way to the interpretation that it may be a signal of new 
physics beyond the SM. Clearly it is still premature to claim that this is a genuine new 
physics effect, mainly due to the large theoretical uncertainties in the non-perturbative 



hadronic sector which affect the SM predictions for e' / e |TJ| . However if one accepts 
the point of view that new CP sources are needed in order to obtain large values for 
e'/e , one may wonder if the minimal supersymmetric extension of the SM (MSSM) 
can help in solving this problem. The answer is no, mainly due to the assumption 
of universal boundary conditions of the soft-breaking terms P, [T5|-|TT|| . Moreover we 
stress that without new flavor structure beyond the usual Yukawa couplings, general 
SUSY models with phases of the soft terms of order 0(1) (but with a vanishing CKM 
phase 5ckm = 0) can not give a sizeable contribution to the CP violating processes 

Iran. 



The impact of the new flavor structure in the non-universality of the A-terms have 
been studied in Refs. • In these works it was emphasized that the non-degenerate 
A-terms can generate the experimentally observed CP violation e and e' / e even with a 
vanishing 5ckm, i-e., fully supersymmetric CP violation in the kaon system is possible 
in a class of models with non-universal A-terms. This effect can be simply understood 



by making use of the mass-insertion approximation : the non-degenerate A-terms 
enhance the gluino contributions to e and e' / e through enhancing the imaginary parts 
of the L-R mass insertions Im(<^ R )i2 and hn{5j tL )x2 @. 

It is well known that the experimental b — > 57 constraints cause a dramatic re- 
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duction of the allowed parameter space in case of universal soft terms pO||2T|. Hence 
one may worry if these constraints would be even more severe in the case of the non- 
degenerate A-terms. However a complete analysis of the b — > 57 constraints has not 
been considered in Refs. |3|-|5|]. In particular in these works it was roughly checked, 
by using the results of Ref. jn5], that the b — > 57 constraints are satisfied, namely 
(^lr)23 < 1-6 x 10~ 2 and (5^)23 < 8-2- Note that these constraints are obtained by 
assuming that the gluino amplitude is the dominant contribution to b — > 57 (dominant 
also with respect to the SM). 

Although the gluino contribution to b —>■ 57 decay is usually very small in the 
universal case, being proportional to the mass insertions (5^)23, we could expect that 
this contribution may be enhanced by the non-universality of the A-terms. However 
the non-universal A-terms could also give large contributions to the chargino amplitude 
through the {S^r)^- Even though the non-degenerate A-terms enhance the gluino 
contribution, one can not a priori expect that it will be the dominant effect. Therefore 
a careful analysis of the b —>■ 57 predictions, including the full SUSY contributions, is 
necessary in this scenario. 

The purpose of this work is to perform a complete analysis of the b —>■ sj constraints 
for the SUSY models with non-universal A-terms studied in Refs. Indeed we 

will show that our results are different from the naive expectations based only on the 
gluino-dominance approximation. 

The paper is organized as follows. In section 2 we present two models with non- 
universal A-terms that have recently been considered for solving the SUSY CP prob- 
lem. These two models are based on weakly coupled heterotic string and type I string 
theories respectively. In section 3 we present formulas for the total branching ratio 
taking into account the QCD corrections and we discuss the different SUSY contribu- 
tions to the b — > S7 decay in these models. Our numerical results for the predictions 
of the b — > S7 branching ratio are presented in section 4. The last section is devoted 
to conclusions. Finally, various formulas are summarized in the appendix. 
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2 String inspired models with non-degenerate 
A— terms 



In this work we consider the class of string inspired model which has been recently 
studied in Refs. In this class of models, the trilinear A-terms of the soft SUSY 

breaking are non-universal. It was shown that this non-universality among the A— 
terms plays an important role on CP violating processes. In particular, it has been 
shown that non-degenerate A-parameters can generate the experimentally observed CP 
violation e and e'/e even with a vanishing <5ckm- 

Here we consider two models for non-degenerate A-terms. The first model (model 
A) is based on weakly coupled heterotic strings, where the dilaton and the moduli 
fields contribute to SUSY breaking |5J. The second model (model B) is based on type 
I string theory where the gauge group SU(3) x U(l)y is originated from the 9 brane 
and the gauge group SU{2) is originated from one of the 5 branes [f|. 

In order to fix the conventions, the following Lagrangian Csb for the soft-breaking 
terms is assumed 

1 1 11 

- C SB = -l',,i.OiOjO h + -{/iB) 13 fafa + -(m 2 )^*^- + -M a \ a \ a + h.c. (1) 

where the fa are the scalar parts of the chiral superfields $j and A a are the gaugino fields. 
In the notation for the trilinear couplings, the A-terms are defined as = Yij^A^k 
(indices not summed) where are the corresponding Yukawa couplings. 

2.1 Model A 

We start with the weakly coupled string-inspired supergravity theory. In this class 
of models, it is assumed that the superpotential of the dilaton (S) and moduli (T) 
fields is generated by some non-perturbative mechanism and the F-terms of S and T 
contribute to the SUSY breaking. Then one can parametrize the F-terms as H 

F s = V3m 3/2 (S + S*)sin9, F T = m 3/2 (T + T*) cos 9. (2) 

Here 717,3/2 is the gravitino mass, rii is the modular weight and tan# corresponds to the 
ratio between the F-terms of S and T. In this framework, the soft scalar masses 
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and the gaugino masses M a are given by || 

m 2 = 7713/2(1 + n>i cos 2 9), (3) 
M a = V3m 3/2 sin9. (4) 

The A u ' d -terms are written as 

(A u ' d ) i:j = -V3m 3/2 sm9 -m 3/2 cos 9(3 + rii + rij +n Hu d ), (5) 

where are the modular weights of the fields that are coupled by this A-term. As 
shown in Eqs.([||)|), the values of the soft SUSY breaking parameters depend on the 
modular weight of the matter fields. These modular weights Tij are negative integers, 
and their 'natural' values (in case of Zn orbifolds) are —1, —2, and —3 ||. If we assign 
77 j = — 1 for the third family and rii = — 2 for the first and second families (we also 
assume that nu 1 = ~ 1 and nn 2 — — 2). Note that with this choice of modular weights 
we have m^ 2 < m 2 Hl which is favored for the electroweak breaking (EW) and all the 
squark mass matrices are equal. Also we find the following texture for the ^-parameter 
matrix at the string scale 



( %u,d -^u,d Vu,d \ 

%u,d •Eu,d Uu,d 
\ Vu,d Vu,d ^u,d / 



(6) 



where 



x u = 77i3/ 2 (— v^sin^ + 3 COS6 1 ), (7) 

Xd = Vu = ^3/2 (-v3 sin 9 + 2 cos 6), (8) 

Vd = z u = m 3/2 (-VSs'm9 + cos9), (9) 

Zd = -V3m 3 / 2 sm9. (10) 

By fixing the value of tan/3 we can determine the values of /x and B from the 
radiatively EW breaking conditions. Then all the SUSY particle spectrum is completely 
determined in terms of 7773/2 and 9. The non-universality of this model is parameterized 
by the angle 9 and the value 9 = n/2 corresponds to the universal limit for the soft 
terms. In order to avoid negative mass squared in the scalar masses we restrict ourselves 
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to the case with cos 2 9 < 1/2. Such restriction on 9 makes the non-universality in 
the whole soft SUSY breaking terms very limited. However, as shown in 0, H], this 
small range of variation for the non-universality is enough to generate sizeable SUSY 
CP violations in K system. We emphasize that choosing modular weights different 
from those assigned above, the allowed range of the soft SUSY breaking terms do not 
essentially change. For instance, if we assign n« = —3 for the first family instead of —2, 
it may appear that the non-universality among the entries of the A-terms is enhanced. 
However in this case the angle 9 is more constrained than before (cos 2 9 < 1/3). We 
have checked that different choices for rii do not significantly affect our results for the 
b — > S7 branching ratio. 



2.2 Model B 

As mentioned in the introduction, this model is based on type I string theory. Like 
model A, this is a good candidate for generating sizeable SUSY CP violations. Recently, 
there has been considerable interest in studying the phenomenological implications of 



this class of models |22| . In type I string theory, non universality in the scalar masses, 
A-terms and gaugino masses can be naturally obtained ||. Type I models contain 
either 9 branes and three types of 5j(z = 1,2,3) branes or 7, branes and 3 branes. 
From the phenomenological point of view there is no difference between these two 
scenarios. Here we consider the same model used in Ref. []5J, where the gauge group 
SU(3)c x U(1)y is associated with 9 brane while SU(2)l is associated with 5i brane. 

If SUSY breaking is analysed, as in model A, in terms of the vevs of the dilaton 
and moduli fields || 

F s = V3m 3/2 (S + S*)sm9, F T ' = m 3 / 2 (T, + T*) e t cos 9 , (11) 

where the angle 9 and the parameter 0j with J2i |©i| — 1> just parametrize the direc- 
tion of the goldstino in the S and Tj fields space . Within this framework, the gaugino 
masses are || 

Mi = M 3 = V3m 3/2 sm9, (12) 
M 2 = V3m 3/2 O lC os9. (13) 
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In this case the quark doublets and the Higgs fields are assigned to the open string 
which spans between the 5i and 9 branes. While the quark singlets correspond to 
the open string which starts and ends on the 9 brane, such open string includes three 
sectors which correspond to the three complex compact dimensions. If we assign the 
quark singlets to different sectors we obtain non-universal A-terms. It turns out that 
in this model the trilinear couplings A u and A d are given by RO] 



A u = A a 



( x y z \ 

x y z 
V x y z J 



(14) 



where 



x = -V3m 3/2 (sm9 + (e 1 -e 3 )cos9), (15) 
y = -y/3m 3/2 (sinO + (Q 1 -e 2 )cos9), (16) 
z = — V3m 3 / 2 sin 6. (17) 

The soft scalar masses for quark-doublets and Higgs fields {ml), and the quark-singlets 
(m^J are given by 



ml = ml /2 ( 1 - |(1 -Ql) cos 2 9), (18) 



= m 2 /2 (l-3e 2 cos 2 fl), (19) 



where i refers to the three families. For 0j = l/y/3 the A-terms and the scalar masses 
are universal while the gaugino masses could be non-universal. The universal gaugino 
masses are obtained at 9 = n/6. 

It is worth mentioning that in these models (A and B) the gaugino masses, the 
A-terms, and the //-term are in general complex. However, by using i?-rotation we 
can make the gaugino masses real and we end up, in addition to the phase of with 
the phases of the A-terms. The phase of fi is severely constrained by the electric 
dipole moment (EDM) of the electron and the neutron WE \, while the phases of the 



A-terms are essentially unconstrained. Thus one can set the phase of fi to zero, as 
done in Ref. |5~Pfl- Moreover it has been shown that the phases of the A-terms can 
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lead to sizeable supersymmetric contribution to CP observables, in particular on the 
direct CP violation e'/e. However, since the total branching ratio b — > 57 decay is a 
CP conserving observable, this should not be very effective in constraining the phases 
of the SUSY soft-breaking terms. For this reason we have made in our analysis the 
simplifying assumption to set to zero all the phases. 

2.3 Yukawa textures 



As emphasized in Refs. in models with non-degenerate A-terms we have to fix 

the Yukawa matrices to completely specify the model. In fact, with universal A-terms 
the textures of the Yukawa matrices at GUT scale affect the physics at EW scale only 
through the quark masses and usual CKM matrix, since the extra parameters contained 
in the Yukawa matrices can be eliminated by unitary fields transformations. This is 
no longer true with non-degenerate A-terms since in the scalar potential the A-terms 
enter through the tensorial product (Y q , )^ = (Y q )ij(A q )ij. Thus the diagonalization 
of Y q can not be done simultaneously with Y A (unlike the universal case). Thus in 
the models with non-universal A-terms, some extra degrees of freedom (in addition 
to the quark masses and CKM matrix) contained in the Yukawa matrices become 
observable. Hence, the analysis of the non-degenerate A-terms could shed some light 
on the favoured Yukawa textures. For instance in Ref. [Q, using a symmetric Yukawa 
texture with a symmetric A-terms, an accidental cancellation between the different 
SUSY contributions to e'/e was found, cancellation which leads to a very small value 
for e'/e. On the contrary, by using asymmetric Yukawa matrices with symmetric A- 
terms, this cancellation does not occur and it is found that e'/e can be easily of the 
order of the KTeV result. 

Here we show two realistic examples of Yukawa matrix textures that have already 
been used in Refs. In the first one we have the following symmetric Yukawa 

matrices 



Y d = y b 



( U 12 ^ V 13 \ [00 v 13 \ 

m, m 



y UOg His y 

12 nib f^b A 

\ V13 U 2 3 1 / 



Y u = y 



q rric q 



(20) 



V v 13 1 J 

where y b,t are the Yukawa couplings of the bottom and top respectively, and V is the 
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CKM matrix. The second example is based on the assumption that the CKM mixing 
matrix originates from the down Yukawa couplings and that the Yukawa matrices are 
hermit ian. 

Y u = — ^— rdiag (m u ,m c ,m t ) , Y d = — — -V f • diag (m d , m s , m b ) ■ V (21) 
v cos p «smp 

Although the analysis of the CP violation is quite sensitive to the specific Yukawa 
matrix, we found that the branching ratio of b — > sj does not essentially depend on 
it. We checked this property by using different Yukawa textures (the two examples 
presented here and others). We will only present through all the paper the results 
concerning the second example just as a representative case. 

Now we present the general expressions for the squark mass matrices in the non- 
universal case in the SCKM basis, in this basis the unitary matrices S R < L and S Dr ' l 
(obtained by a superfield rotation) are chosen to diagonalize the up- and down- Yukawa 
couplings Y u,d 

vsm/3 / vu ,dxTr,t _ _ vcosB d T f 

m u - ^ b u R [i ) o Ul , m D - — &D R {* ) Dl 

where T stands for the transpose and rnu,D are the diagonal up- and down-quark 
mass matrices respectively. In this basis the up and down squark mass matrices at low 
energy (respectively M| and Mj) are given by p3] 




where for the up-sector 

2 

(Ml) LL = S UL M 2 Q Sl L + ml + ^(3- 4sin 2 W ) cos 2/3, 

( m *)rr = s UR {Ml c fs\ jR + m l + 2 -^^e wC o^ 

KL = (Ml)^ = »m uC otp+ V ^S UL Y u A *Sl R , (24) 
and for the down-sector 

2 

(Mj) LL = S Dl MIS ] Dl + ml- ^(3 - 2 sin 2 M cos 2/5, 
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(Mj) RR = S DR {Ml) T S ] DR + ml + ^sm 2 6 wC os2P 

v cos (3 
72 
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= {Mj)l L = »m D UnP + V ^S DL YtSh R (25) 



where M| and M? c ^ c are the soft-breaking (3 x 3) mass matrices for the squark doublet 
and singlets respectively. Our convention for the sign of the /i-term is opposite to the 
same one of Ref. ||23|| . Note that the matrices Sv,d, unlike in the universal case, can 
not be re-absorbed in the definition of diagonal Yukawa couplings. 



3 The b —> 57 decay in SUSY models 

In this section we analyse the b — > S7 decay in SUSY models with non-universal A— 
terms. As pointed out previously, these models are particularly interesting because 
they can give sizeable contributions to the CP violating processes through their large 
contributions to the mass insertions {S^p)ij or (6^)ij. For this reason one may expect 
that large effects can be also induced in the processes mediated by the dipole-magnetic 
operators, such as the rare decay b — > 57 . Indeed, if the S^ R are large enough, 
then the SUSY contributions to these operators are enhanced since the typical chiral 
suppression is removed by the insertion of the internal gaugino mass, thus allowing for 
a competition with the chiral-suppressed SM amplitude. 

Let us start with the experimental results. The most recent result reported by 
CLEO collaboration for the total (inclusive) B meson branching ratio B — > X s 7 is p4 



BR(5 -> X s7 ) = (3.15 ± 0.35 ± 0.32 ± 0.26) x 10" 4 (26) 

where the first error is statistical, the second systematic, and the third one accounts 
for model dependence. From this result the following bounds (each of them at 95% 
C.L.) are obtained 

2.0 x 10~ 4 < BR(5 -> X sl ) < 4.5 x 10" 4 . (27) 
In addition the ALEPH collaboration at LEP reported a compatible measurement of 



the corresponding branching ratio for b hadrons at the Z resonance |25 . 

The starting point for the theoretical study of b — > 57 decay is given by the effective 
Hamiltonian 

AC 8 

Heff = y=- V* 2 V 33 £ CMQM (28) 

i=l 
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where the complete basis of operators in the SM can be found in Ref. ||26|| . Recently 
the main theoretical uncertainties present in the previous leading order (LO) SM cal- 
culations have been reduced by including the NLO corrections to the b — > sj decay, 
through the calculation of the three-loop anomalous dimension matrix of the effective 
theory (EBl. The relevant SUSY contributions to the effective Hamiltonian in Eq. 
affect only the Qj and Q$ operators, the expression for these operators are given (in 
the usual notation) by 

e 



Qr 
Qs 



16vr 2 

9s 



m b {s L a^b R ) Ffu, , 
m b (s L a^T a b R ) G° 



(29) 



The Wilson coefficients Cj(/i) are evaluated at the renormalization scale //& ~ 0(m b ) 
by including the NLO corrections |26[. They can be formally decomposed as follows 

(30) 



An 



where Cf^ and C\ ^ stand for the LO and NLO order respectively. Finally the branching 
ratio BK(B — > X s j), conventionally normalized to the semileptonic branching ratio 
BR exp (B -> X c ev) = (10.4 ± 0.4)% 0, is given by |B 



BR NLO (£ -> X s7 ) 



x 



BR exp (B -> X c ez/ 



1^33 


2 6a e 




^23 


2 


7T£(z)A;(z) 



1 3 



a s (m b ) 



71 



npj j 



(31) 



with 



«s(^6 
47T 



(0), 



(M 6 )log<5(7+21og<S)/37r 



i=l 



l)|CfW + 



r i(z) +7 4 '? lo . 



a s (yWfc 



7T 



E Q^m^mms) 
i<j=i 



where z = m 2 c /m1. The expressions for Cl°\ C\ , and the anomalous dimension matrix 
7, together with the functions g(z), k(z), r^z) and fij(5), can be found in Ref. 



The term 5 np (of order a few percent) includes the non-perturbative l/m^ [28| and 
l/m c P9j corrections. From the formula above we obtain the theoretical result for 
BR(i? — > Xg'-f) in the SM which is given by 



BR NLO (5 -> X sl ) = (3.29 ± 0.33) x 10" 



(32) 
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where the main theoretical uncertainty comes from uncertainties in the SM input pa- 
rameters, namely m t , a s (Mz), ot emi m c /mb, m&jV^, and the small residual scale de- 
pendence. | The central value in Eq. (|32| ) corresponds to the following central values for 
the SM parameters m v t oXe ~ mf(m z ) ~ 174 GeV, mf = 4.8 GeV, mP c olc = 1.3 GeV, 
fi b = m b , a s (mz) = 0.118, a f T 1 (m^) = 128, sin 2 ^ = 0.23 and a photon energy res- 
olution corresponding to 5 = 0.9 is assumed. Note that in Eq.(^Tj) the (small) l/m c 
corrections have not been included. 

The SUSY contributions to the Wilson coefficients are obtained by calculating 
the b — > sj and b — > sg amplitudes at EW scale respectively. The LO contributions to 
these amplitudes are given by the 1-loop magnetic-dipole and chromomagnetic dipole 
penguin diagrams respectively, mediated by charged Higgs boson, chargino, gluino, and 
neutralino exchanges. The corresponding results for these amplitudes can be found in 



Ref. [21]. It is known that the charged Higgs contribution always interferes with the 
SM contribution JZI], The chargino contribution could give rise to a substantial 
destructive interference with SM and charged Higgs amplitudes, depending on the sign 
of /i, the value of tan/3 , and the mass difference between the stop masses [PU|,PT|]. 

We point out that the SUSY models with non-universal A-terms may induce non- 
negligible contributions to the dipole operators Qy^ which have opposite chirality with 
respect to Qt^. It is worth mentioning that these operators are also induced in the 
SM and in the MSSM with supergravity scenario, but their contributions are negligible 
being suppressed by terms of order 0(ra s /m&). In particular in MSSM, due to the 
universality of the A-terms, the gluino and chargino contributions to Qj^ turn out to 
be of order 0(m s /m&). This argument does not hold in the models with non-universal 
A-terms and in particular in our case. It can be simply understood by using the mass 
insertion method For instance, the gluino contributions to Q-j and Qi operators 
are proportional to (<^) 23 ~ {S Dl Y^S ] DR ) 2 ^/m] and (5 d RL ) 23 ~ (S DR Y^S ] DL ) 2Z /m] 
respectively. Since the A D matrix is symmetric in model A and Af- ~ A® in model B, 
then (Sf jR )23 — (8rl) 23- Then in our case we should consistently take into account the 



§ Recently in Ref. |30| ] the current method of extracting the inclusive rate for b — > sj from the 
currently published CLEO data has been criticized arguing that the theoretical uncertainties have so 
far been underestimated and only a precise measurement of the photon spectrum would be help in 
reducing these uncertainties. 
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SUSY contributions to Q 7 in b — > s-y . Analogous considerations hold for the operator 

Qs- 

By taking into account the above considerations regarding the operators (37,8, the 
new physics effects in b — > 57 can be parametrized in a model independent way by 
introducing the so called R 7 ^ and i?7,8 parameters defined at EW scale as 

/WO) r (0)SM\ ~( ) 

^7,8 - (0)SM ' ^7,8 - ( )SM ' W 

^7,8 ^7,8 

where C 7) 8 include the total contribution while Cf$f contains only the SM ones. Note 
that in CV,8, which are the corresponding Wilson coefficients for Q 7 $ respectively, we 
have set to zero the SM contribution. In Ref. |21j only the expressions for the R 7j $ 
are given, for completeness we report the corresponding expressions for R 7 8 in the 
appendix. 

Inserting these definitions into the BK(B — > X s 7) formula in Eq.(|3lD yields a 
general parametrization of the branching ratio in terms of the new physics contributions 

HQ 



BR(5^X S 7) = (3.29 ± 0.33) x 10 -4 (l + 0.622i2 7 + 0.090{R 2 7 + R 2 7 ) 

+ 0.066^8 + 0.019(^8 + RyRs) + 0.002(R 2 8 + R 2 8 )) , (34) 



where the overall SM uncertainty has been factorized outside. We have checked ex- 
plicitly that the result in Eq.(|34]) is in agreement with the corresponding one used in 
Ref. |32| . Recently the leading EW corrections in the SM have been included in the C 7 $ 



coefficients ||33[| . However, we have not included them since they should affect the re- 
sults in a few percent which is within the theoretical uncertainties present in the SUSY 
sector. In particular, in order to be consistent with the NLO calculations, one should 
also include the corresponding two-loop QCD corrections to the SUSY amplitudes, 
namely C 7 l \. In this respect recent works [[34J have been done in this direction. They 



calculated the two-loop a s corrections to the chargino and gluino amplitudes. However 
these corrections can be consistently applied only in a restricted SUSY scenario with 
low tan p and large gluino masses, since on the contrary other two-loop diagrams with 



^Note that the SM central value for BK(B — > X s "f) in Ref. |3l]] slightly differs from the result in 
Eq.(|3 7 ^), since in Ref. |H| the non-perturbative A/m c corrections have been included. 
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no QCD interactions at all could become relevant. In our work we are going to explore 
SUSY scenarios in the full range of tan (3 (2 < tan/? < m t /mb) and gluino mass, so we 
do not include these corrections since they would not change our main conclusions. 



4 Numerical results and discussions 

In this section we present our results for the total branching ratio BK(B — > Xg'y) 
as a function of the fundamental parameters of the soft-breaking sector of models A 
and B. We start our discussion by analysing the constraints on these models set by 
the condition of vacuum stability and the experimental bounds on the SUSY particle 
spectrum. 

We calculate the low energy SUSY spectrum by running the soft-breaking terms and 
the other SUSY parameters (by means of the renormalization group equations (RGE) 
of MSSM generalized to the non-universal soft-breaking terms [EH]]) from the GUT scale 
(~2x 10 16 GeV) to the EW scale (~ Mz), by using the boundary conditions given in 
section 2. For fixed tan /3 and sign of fi, we restrict the parameter space by imposing 
the present experimental bounds on the SUSY spectra PTfl . Since the spectrum of 



these models is given in terms of few parameters (m 3 ^, 9 and/or 0j) we find that by 
requiring the lightest chargino mass to be m x ± > 90 GeV implies that all the other 
SUSY particle and the lightest Higgs masses are above their experimental bounds. 

In order to avoid vacuum instabilities and color-charge breaking, we require that 
all the square scalar masses in Eqs.(||), ( |T9"D should be positive. All these constraints 
set strong restrictions on the parameter space of both models. As pointed out in 
section 2, in model A this leads to sin 9 > 1/V2 ~ 0.7. The bounds on m 3 / 2 from the 
lightest chargino mass are as follows: for low tan/3 (tan/3 = 2) m 3 / 2 > 80 (60) GeV for 
s'm9 ~ 1/V2 (1). For large tan/? these bounds are increased by 20 GeV in both cases. 

In model B the non-universality is parameterized by the angle 9 and the 0j's. The 
values 0j = 1/ -\/3 give universal A terms and scalar masses, the gaugino masses are 
universal only at 9 = n/6. To avoid negative scalar masses in this model one needs to 
impose the constraints 



cos 2 9 0^ < 1/3, (i = 1, 2, 3) , 
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(l-6?)cos 2 # < 2/3. 



(35) 



In the limit of universal A terms these two constraints are satisfied for any value of 9. 

Further constraints on the parameter space of model B are obtained from the gaug- 
ino sector. As shown in Eq. (p~3|) , for a non very large 1713/2, the limit of 9 ~ 71/ 2 can 
not be reached since M2 would approach zero and the mass of the lightest chargino 
would be too small. Moreover, the lower bound on the gluino mass and the condition 
of having the EW breaking at the correct scale, lead also to a lower bound on the angle 
9. In the case of universal A-terms we find that 9 > 0.1, while in the non-universal 
case Eqs. (|35D set more severe constraints on 9. For example, for 0i ~ 1 we obtain that 
0.9 < 9 < 7r/2 and for 9 close to 7r/2 the gravitino mass m.3/2 should be quite heavy (of 
order TeV) to make the lightest chargino higher than the experimental bound. We find 
that these constraints together strongly reduce the allowed parameter space of model 
B. 

We have checked that, in both models, the B — B mixing measurements do not set 
further constraints on the allowed ranges of the parameter space. This mixing, being 
a AB = 2 process, is proportional to (<5l#)i3 x (5^£»)i3 where A, B,C,D = (L, R). We 
have found that, in our case, the values of these mass insertions satisfy the constraints 
given in Ref . . 



Our results for the partial SUSY amplitude contributions are presented in Figs. [|T]-U . 
The total branching ratio BK(B — > Xg'j) is shown in Figs. and for model A and 
B respectively. In Figs.|l|-|3[ we show, for model A, the individual SUSY contributions 
to the i?7 variable, (see Eq.(^) for its definition) versus sin#. We see that the chargino 
and charged Higgs contributions give the dominant effect in all the range of 9, while the 
gluino is sub-dominant, such as in the universal case. For model B, as expected from 
its tightly constrained parameter space, we have found that the results of the separate 
amplitude contributions do not differ from the corresponding ones in the universal 
scenario, and will not be presented here. 

We can understand this behavior by using the mass insertion method. The gluino 
amplitude gets two leading contributions: one is proportional to the single-mass- 
insertion (6f jR )23 and the another one to the double-mass-insertion, namely (5^)22(^^)23- 
In the low and intermediate tan f3 regions these two mass insertions are comparable, so 

14 



a possible destructive or constructive interference between them may appear depending 
on the sign of fi. In the large tan /3 region the double mass insertion becomes dominant, 
since (5^ H ) 22 is proportional to //tan/5 and the amplitude (normalized to the SM one) 

-P- tan/3 (m 2 w M~ g fa) (36) 

where rriq is the average squark mass in the down sector and M~ _ is the off diagonal 
element of the down-squark mass matrix defined in Eqs. (|23]j25| ). In this way, the 
sensitivity of the gluino amplitude to the sign of fi for large tan (3 can be understood. 
The chargino amplitude, like in the gluino case, is enhanced by tan /3 and the dominant 
contributions to it are the (Higgsino) A~ h _ and (Wino-Higgsino) A^- amplitudes, given 
by 

Ar Mf r 

— oc tan/3 U m t ) * Lt \ , (37) 

A SM m t L m t R 

A - - 1 M~ - 

_pz. oc tan/?(/xM 2 ) -ft. (38) 

V32 m t L m t R 

We see that these contributions are large with respect to the gluino one, such as in the 
universal case of MSSM, mainly because the mass insertions get a large enhancement 
of the light-stop mass m| in the denominator, unlike in the gluino case where the 
l/m| suppression is effective. 

Regarding the contribution of the Qj^ operators to the total branching ratio, we 
checked that their effect is negligible, almost an order of magnitude smaller than the 
total contribution to Qj t s- This can be explained by observing that the dominant effect 
to these operators comes from the gluino amplitude which is much smaller than the 
chargino or charged Higgs one. 

In Figs. we plot the results for the branching ratio BR(Z? — ► -XVy), in model 
A, versus sin# for different values of tan/3 , the sign of fx and for two representative 
values of gravitino mass 777.3/2 , namely 7773/2 = 150, 300 GeV. The main message arising 
from these results is that the sensitivity of BK(B — > X s 7) respect to sin# increases 
with tan/3 . In particular for the low tan/3 region the b —> sj result does not differ 
significantly from the universal case. In the large tan/3 region, tan/3 = 15 — 40, the 
CLEO measurement of b — >• 57 set severe constraints on the angle 9 for low gravitino 
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masses. For /j < almost the whole range of parameter space is excluded as shown in 
Fig- §§■ 

Comparing Top and Bottom plots in Figs.@||, we see that that for positive (neg- 
ative) sign of fi the branching ratio BR,(B — > X s 7) decreases (increase) when the 
departure from the universality increases (sin# — > 1/^/2). Clearly, due to decoupling 
effects, the deviations from universality tend to be reduced for large gravitino masses, 
as can be seen by comparing the plots at 7713/2 = 150 GeV with the corresponding ones 
at m 3 /2 = 300 GeV. 

Now we discuss results in model B. In Fig.|3j we plot the branching ratio BK(B — > 
Xs'y) versus tan /3 for three different values of 0i, 02 (see the figure caption) which are 
representative examples for universal and highly non-universal cases. From these fig- 
ures it is clear that BR(I? — > X s ^) is not very sensitive to the values of Oj's parameters, 
even at very large tan j3 , unlike model A. The constraints from CLEO measurement 
are almost the same in the universal and non-universal cases. For fi > the branching 
ratio is constrained from the lower bound of CLEO only at very large tan (3 , while for 
/j < the branching ratio is almost excluded except at low tan /3 . 

5 Conclusions 

The recent CP violation measurements of e'/e indicate a large deviation from the SM 
predictions which might be interpreted like a signal of new CP violation sources beyond 
the SM. It is unlikely that the minimal supersymmetric standard model with universal 
boundary conditions at GUT scale can explain these large enhancements in the direct 
CP violations. On the contrary, non-minimal SUSY scenarios with non-universal A— 
terms, derived from some string inspired models, have been found effective in explaining 
large values for e'/e while keeping the electric dipole moments below the experimental 
bounds. 

In this paper we have considered two models based on weakly coupled heterotic 
string (model A) and type I string theories (model B). In this framework we carefully 
analysed the constraints set by the b — > 37 decay on these two models by taking into 
account the relevant set of SUSY diagrams. In the calculation of the total branching 
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ratio we take into account the NLO QCD corrections to the SM. We found that in 
the model based on the weakly coupled heterotic string, the b — > 57 branching ratio 
is more sensitive to the non-universality at large tan/? and that the dominant SUSY 
contribution comes from the chargino amplitude for any value of tan/3 . For type I 
string-derived model, we found that the sensitivity of the b — *■ sj branching ratio to 
the non-universality parameters 9 and Gj is quite weak. The main reason for this 
weakness is because in this model the allowed ranges for these parameters are strongly 
constrained by the vacuum stability bounds and the experimental limits on the lightest 
chargino mass. 

We conclude that the recent CLEO measurements on the total inclusive B meson 
branching ratio BR(i? — > X s j) do not set severe constraints on the non-universality of 
these models. Moreover the constraints set on tan /3 and gravitino mass are almost the 
same as in the universal case. In this respect we have found that the parameter regions 
which are important for generating sizeable contributions to e'/e in particular 

the low tan (3 regions, are not excluded by b — > 57 decay. 
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Appendix 

Here we give the expressions for the dominant SUSY contributions to Rj^ defined in 
Eq. (|5B|) , namely the chargino and gluino ones, in the approximation 0(m s /m&) = 

R7.8 — + -^7,8 



3V 3 * 2 V 33 x tw F 7 (x t w) 
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with 



x 



1=1 k=l 



_"Xi 

2 6 



m 



m 



3X^32 V 33 xtwFi(xtw) J^kTi 
16a s 6 



* m 
fc2 m& 



£*™,( rD 'L( rD ")LS F4( ^ 



27a w V3 2 V 3 3X W F 7 (x tT y) ^ 
2« 5 



9«m/ P32 ^33x^^1 (x tw ) 



X 



fe=l 



A' 



-v n (r^ y) + 1 , - v J2 (r^ v) 

-— ^ -U /2 (r^ V M D ) 

\/2m w cos /? v 7 fci 



(39) 



(40) 



where F 7 (x) = jFi(z) + F 2 (z), U and V are the 2x2 diagonalization 

matrices for the chargino mass matrix defined as in Ref. ||21]|, and the 6x6 matri- 



ces 



T {u,D) = 



r (D L ,U L ) r (D R ,U R ) 
1 6x3 5 1 6x3 



respectively diagonalize the Up- and Down-squark 
mass matrices in Eqs. (p3| , |25| ). The expressions for the functions Fi(x) can be found in 
Ref. El. 
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Figure 1: (Top) Rj for chargino contribution versus sin# in model A, for fi > 0, 
m 3/2 = 150 GeV, and for tan/3 =2, 15, 40. (Bottom) The same as before, but for 
H < 0. 
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2: (Top) R-j for the gluino contribution versus sin# in model A, for fi > 0, 
: 150 GeV, and for tan/3 =2, 15, 40. (Bottom) The same as before, but for 
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Figure 3: R 7 for charged Higgs contribution versus sin 6* in model A, m 3 / 2 = 150 GeV, 
and for tan/5 =2, 15, 40. 
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Figure 4: (Top) The BR(B — > X s ^) versus sin# in model A, for fj, > 0, m 3/2 = 150 
GeV and tan/5 =2, 15, 40. (Bottom) The same as before, but for /i < 0. 
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Figure 5: (Top) The BR(B — > X s ^) versus sin6» in model A, for fx > 0, m 3 / 2 = 300 
GeV and tan/3 =2, 15, 40. (Bottom) The same as before, but for /i < 0. 



26 



5 
4.5 



o 








_><, 


3.5 - 






03 


3 - 


X 




A 






2.5 - 


CO 




DC 




CD 






1.5 - 



i i i i i i i 

CLEO 




SM 




CLEO 


i i i i i i i 



10 

9 

^f~~ 8 

1 

o 

7 

><^ 

'P 6 

03 
X 

A 



CO 

DC 
CQ 



10 



15 



20 25 

tan(P) 



30 



35 



40 



I I I 


1 1 1 1 ^ 




CLEO " 




SM 




CLEO 


I I 1 1 I I 1 



10 



15 



20 25 

tan(P) 



30 



35 



40 



Figure 6: (Top) The branching ratio BR(B — > X s ^) versus tan (3 in model B, for /i > 0, 
m 3/2 = 150 GeV, and for some values of (6i, 6 2 ) = (1/a/3, l/>/3), (0.9,0.2), (0.6,0.2), 
corresponding to the continuos, dashed, and dot-dashed lines respectively. (Bottom) 
The same as before, but for ji < 0. 
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